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ABSTRACT 


A  summary  of  microwave  gas  breakdown  theory  is  presented  with 
special  emphasis  on  nanosecond  pulse  breakdown.  The  various  approaches 
previously  presented  in  the  literature  are  compared  and  shown  to  be  identical 
as  far  as  practical  results.  It  is  suggested  that  in  general,  for  pressure-time 
products  greater  than  1  torr-nanosecond,  measurements  made  under  cw 
conditions  will  be  satisfactory,  A  program  of  future  investigation  is  outlined 
to  test  this  suggestion  and  to  gather  the  needed  data  on  pulse  breakdown. 
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NANOSECOND  PULSE  BREAKDOWN  STUDY 


I.  INTRODUCTION 

This  report  summarizes  a  study  of  breakdown  initiated  by  intense 
short -duration  microwave  pulses.  We  will  present  the  results  of  our  analysis 
concerning  prior  theoretical  and  experimental  approaches  and  then  give  a 
critic  of  these  methods.  The  importance  of  individual  processes  in  a  gas 
that  effect  breakdown  will  be  discussed.  A  future  program  is  also  outlined 
that  incorporates  a  combined  empirical  and  theoretical  attack  on  the  problem. 

The  problem  of  breakdown  for  nanosecond  pulses  arises  out  of  the  need 
for  higher  peak  pulse  power  in  shorter  pulses  to  increase  the  range  and  range 
resolution  of  present  radar  systems.  Previous  studies  of  pulsed  microwave 
breakdown  have  been  directed  primarily  towards  the  problem  of  satellite 
communication.  The  radar  detection  problem  is  somewhat  different  in  that 
the  equipment  is  all  ground  based.  Thus,  because  of  the  well  known  variation 
of  breakdown  with  pressure,  a  breakdown  minimum  will  occur  at  an  altitude 
of  about  100,  000  ft.  That  is,  unless  some  special  focusing  scheme  is  used 
we  may  expect  breakdown  to  occur  in  the  pressure  range  of  1  atmosphere  to 
about  1  torr.  To  allow  for  pressurized  equipment  in  our  analysis  we  will 
consider  pressures  from  1  torr  to  6  atm.  This  is  a  convenient  division  since 
it  turns  out  that  in  the  low  pressure  region  a  different  mechanism  must  come 
into  play.  We  will  consider  frequencies  in  the  gigacycle  range,  although  as 
will  be  shown  this  is  not  critical.  The  analysis  will  apply  to  pulse  lengths 
of  1-1000  nanoseconds. 

The  general  plan  is  to  first  discuss  the  basic  processes  of  breakdown 
including  the  generation  and  loss  of  electrons,  criteria  for  breakdown,  and  the 
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interaction  of  the  e.  m.  wave  with  the  ionized  gas.  We  shall  discuss  the 
processes  and  their  relative  importance  for  short  pulse  breakdown.  After 
reduction  of  the  loss  processes  by  physical  arguments  to  attachment  we  shall 
derive  the  equations  that  govern  breakdown.  Microwave  testing  methods  and 
future  experiments  are  also  discussed. 
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II.  BASIC  PROCESSES  OF  BREAKDOWN 


1)  General  Considerations 

We  are  considering  the  problem  of  breakdown  in  free  air,  that  is,  in 
a  medium  not  bounded  by  a  container.  This  greatly  affects  the  type  of  analysis 
we  shall  apply  to  the  problem.  In  addition  we  are  interested  in  breakdown 
because  it  affects  the  propagation  of  the  microwave  energy.  The  problem  we 
are  examining  is  a  general  one  of  the  interaction  of  an  e.  m.  field  and  an 
ionized  gas.  This  problem  arises  in  communicating  with  a  re-entering  vehicle 
surrounded  by  a  thermally  ionized  gas,  in  plasma  diagnostics,  and  in  the 
propagation  of  signals  strong  enough  to  cause  ionization.  The  general  character 
of  the  problem  should  be  kept  in  mind  both  in  defining  breakdown  and  in  searching 
the  literature. 

In  general  use  the  term  breakdown  indicates  failure  of  a  normally 
insulating  substance  so  that  it  becomes  conducting.  In  a  gas  this  means  the 
presence  of  free  electrons  obtained  by  ionization  of  the  atoms.  In  d.  c.  or 
Townsend  measurements  one  may  use  the  free  electron  density  as  a  direct 
indication  of  breakdown,  that  is  one  may  observe  the  temporal  growth  of 
the  current  density.  For  microwave  frequencies  this  direct  measurement  is 
neither  convenient,  nor  simple.  The  conductivity,  which  is  directly  related 
to  the  free  electron  density  will  strongly  affect  the  propagation  of  electro¬ 
magnetic  waves  in  the  gas.  We  may  therefore  use  the  wave  causing  breakdown 
as  a  direct  measuring  device  to  determine  the  electron  concentration.  We 
must  define  an  ionization  level  that  we  may  use  as  a  breakdown  criterion. 
Sometimes  an  arbitrary  electron  density  or  an  arbitrary  ratio  of  final  to 
initial  densities  is  used.  We  will  use  a  density  that  is  directly  connected  to 
the  propagation  properties  of  the  wave.  This  is  the  "resonant  plasma  density", 
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which  may  be  defined  as  the  electron  density  n  at  which  the  plasma  angular 


frequency^Jp  is  equal  to  the  applied  angular  frequencyCO  . 
quency  is  (Ref*  1) 


The  plasma  fre- 


2  =  4Tfe2' 
U)  p  m 


n  , 


(1) 


And  we  may  solve  for  the  breakdown  density  n^  in  terms  of  the  signal 
frequency. 


n,  =  1.  24  x  10"8  f2  , 

D 


(2) 


where  e  and  m  are  the  electron  charge  and  mass  respectively. 

We  shall  show  how  this  relates  to  the  propagation  of  microwaves.  As  a 
simple  approximation  consider  the  case  of  a  wave  normally  incident  on  an 
abrupt  air-plasma  interface.  The  reflection  coefficient  R  for  this  case  in 
terms  of  a  generalized  index  of  refraction y\  is 


<3> 

For  a  homogeneous  isotropic  plasma  we  have  (Ref.  2) 

„2  = , .  JsLeVl  ,4, 

where  the  collision  frequency  of  the  electrons.  Figure  1  shows  the 

magnitude  of  the  reflection  coefficient  R  as  a  function  of  fp/f.  We  see  that 
in  the  collisionless  case  there  is  an  abrupt  reflection  at  f  =  f^.  The  other 
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cases  are  not,  so  sharply  affected  and  appreciable  transmission  occurs  for 

f  >  f.  However,  we  have  thus  far  not  included  attenuation  which  also  affects 
P 

the  propagation.  The  breakdown  criteria  given  above  has  also  been  used  by 
Margenau  and  others. 

We  will  now  study  the  processes  that  govern  the  growth  of  the  free 
electron  density  in  a  gas.  In  a  high  frequency  field  an  electron  will  oscillate 
back  and  forth,  and  because  of  inertial  effects  will  be  out  of  phase  with  the 
field.  The  electron  will  not  gain  energy  during  these  oscillations,  unless  it 
happens  to  collide  with  a  gas  atom.  When  an  elastic  collision  occurs  the 
electron,  because  of  its  very  much  smaller  mass,  will  keep  most  of  its  energy, 
but  will  be  randomly  scattered  with  a  different  phase  angle.  It  will  continue 
to  oscillate  back  and  forth  until  it  has  another  collision.  If  the  number  of 
collisions  per  cycle  is  high  the  electron  can  gain  appreciable  energy  from  the 
field.  Since  the  motion  is  randomized  the  electron  gas  is  in  effect  heated  up. 

The  excursion  of  the  electrons  in  a  high  frequency  field  is  rather  small 
so  that  the  effects  are  confined  to  the  volume  of  the  gas,  and  electrode  effects 
will  be  small.  On  the  other  hand,  in  a  d.  c.  discharge  the  electrons  are  swept 
out  of  the  gas  and  breakdown  is  controlled  by  secondary  effects  at  the  electrodes 
and  walls.  (Ref.  3)  This  may  also  happen  for  lower  frequency  r.  f.  discharges. 
(Ref.  4)  In  high  pressure  discharges  there  may  also  be  a  transport  of  photons 
to  the  electrodes  causing  photoelectric  emission.  However,  we  shall  concern 
ourselves  here  only  with  volume  effects.  We  now  write  an  equation  for  the 
creation  and  destruction  of  free  electrons.  Electrons  will  be  created  only  by 
ionization  in  the  volume.  We  may  break  this  into  two  terms:  a  constant  due  to 
cosmic  rays  or  other  natural  sources  plus  a  part  due  to  collisional  ionization. 
Electrons  may  be  lost  by  diffusion,  by  attachment  to  form  stable  negative  ions, 
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and  by  recombination  with  the  positive  ion*.  The  general  equation  governing 
the  electron  density  is  then 

-fr  =^in  '  'Kn -o(n2  -VT+s  (5) 

where  ^/.  is  the  number  of  ionizing  collisions  per  electron  per  second,  "JS 
i  a 

is  the  attachment  coefficient,  (X  is  the  recombination  coefficient,  which  depends 
on  the  product  of  the  electron  and  positive  ion  densities,  which  we  assume  are 
the  same,  although  this  may  not  be  correct  if  the  negative  ion  density  is 
appreciable.  ^7*  is  the  divergence  of  the  diffusion  current  density  and  S  is 
a  general  source  term.  For  example,  S  may  represent  the  source  term  due 
to  x-rays  introduced  to  eliminate  statistical  variation  in  the  measurement 
during  an  experiment. 

It  can  easily  be  seen  that  recombination  is  negligible  for  small  ion 

densities,  especially  if  attachment  is  important.  Loeb  (Ref.  5)  gives 

°(air  ^  10  ^  cm^/ion  sec.  This  leads  to  a  loss  rate  of  2  ions  per  cm^ 

6  3 

per  microsecond  if  n  is  10  /cm  .  Attachment  under  the  same  conditions 
3  3 

will  be  about  10  p  per  cm  per  microsecond,  where  p  is  the  pressure  in  torr. 
For  these  low  densities  the  diffusion  rate  will  also  always  be  much  greater 
than  the  recombination  rate.  In  the  afterglow  when  the  electronic  field  is 
removed  recombination  may  be  important;  however,  for  electron  attaching 
gases,  attachment  may  always  be  dominent  because  attachment  cross-sections 
are  usually  high  for  low  energy  electrons,  and  in  the  afterglow  we  expect 
thermal  electrons.  In  C.  W.  breakdown  Herlin  and  Brown  (Ref.  6)  have  shown 
that  diffusion  is  the  dominant  loss  mechanism,  but  in  our  case  especially  at 
higher  pressures  we  shall  show  that  diffusion  is  negligible. 
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2)  Diffusion 


In  this  section  we  shall  discuss  the  processes  of  diffusion  and  the 
importance  of  diffusion  as  a  loss  mechanism.  It  will  be  shown  that  in  the 
pressure  time  regime  under  consideration  diffusion  is  a  negligible  loss 
mechanism. 

As  discussed  in  Appendix  I  t.he  particle  current  density 

r  =  -vu>»>  .  (6) 

where  D  is  the  diffusion  coefficient  and  is  a  measure  of  the  mean-squared' 

distance  a  particle  diffuses  per  unit  time.  This  may  be  substituted  into 

equation  (5),  Where  ^  is  considered  negligible,  and  v-  w-y.-  This  reduces  to 

1  a 

=  -pn  +  D y2  n  +  S  (7) 

if  we  assume  that  D  is  not  a  function  of  r.  General  solutions  to  this  equation 

can  be  found,  but  they  are  not  to  useful  because  they  depend  on  the  geometry, 

and  except  in  the  most  elementary  cases  they  are  extremely  difficult  to  obtain. 

We  can  obtain  all  the  information  we  require  by  realizing  that  the  solution 

will  involve  a  characteristic  diffusion  length  (/\)  which  will  be  of  the  order  of 

the  actual  physical  dimensions  of  the  system  (i.  e.  —  for  infinite  parallel 

plates  with  spacing  L).  Therefore  we  may  calculate  a  characteristic  diffusion 
2 

time  (/^  /D),  and  assume  that  diffusion  is  negligible  if  the  pulse  length  is 
much  smaller  than  the  diffusion  time.  We  have  shown  in  the  appendix 
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that 


D  = 


2 

Jm 


(8) 


where  u  is  the  energy  of  the  electrons,  and^  is  the  collision  frequency. 

In  general  is  a  function  of  the  velocity.  The  variation  of  ~JJ ^  with  velocity 
will  depend  on  the  interaction  potentials.  The  exact  solution  of  this  problem 
is  extremely  difficult,  however  there  are  two  approximations  often  made  that 
lead  to  simple  results.  These  are,  hard  sphere  collisions  and  an  inverse 
fifth  power  repulsion  law,  which  lead  to  constant  collision  frequency  (l^c)  and 
constant  mean  free  path  (^ )  respectively.  While  it  is  not  expected  that  either 
of  these  assumptions  accurately  describes  the  interaction  in  a  gas  they  have 
provided  reasonable  solutions  for  many  of  the  common  phenomena.  Either  of 
these  assumptions  will  yield  a  similar  result  for  D. 

Assume  that  vc  is  independent  of  velocity  then  we  find  for  air  that 


D  =  2.  2  x  10 5 


(9) 


9 

where  (  u^  is  in  electron  volts  and  we  have  used  =  5  x  10  p  as  given 
by  Brown.  (Ref.  8)  We  may  also  write 


D  = 


(10) 


The  brackets  (  )  denote  an  average  over  the  velocity  distribution  function 
in  this  case  J{P), 
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which  after  substitution  oi  X  =  v/  ^  and  assuming  J(  is  constant  gives 

d  =  t^(v>  <n> 

<V> 

The  mean  free  path  is  just  £  .  If  we  use  our  value  for  y  we  may  compare 

c  < 


these  two  diffusion  coefficients 


D=  <v2)  _  (v2) 


for  a  Maxwellian  distribution  this  ratio  is  1.2.  We  seek  a  value  for  ^  u^  . 

This  depends  on  the  gas  and  the  fields  but  it  certainly  is  not  greater  then  the 
ionization  potential  of  say  10  ev  for  air.  So  that  D  =  and  for  pf\  >  1 

we  have  diffusion  times  in  the  order  of  a  microsecond.  We  may  therefore 
neglect  diffusion  in  our  theory  of  nanosecond  pulse  breakdown.  If  one  wants 
to  consider  extremely  low  pressure  or  perhaps  an  experiment  where  spacings 
are  made  very  small,  then  diffusion  must  be  taken  into  account.  For  repetitive 
pulses  diffusion  as  well  as  recombination  is  important  in  the  afterglow  between 
pulses. 

These  approximate  diffusion  coefficients  depend  on  the  average  electron 
energy,  which  is  a  function  of  E/p  .  In  a  gas  discharge  the  electron  distribution 
is  altered  from  the  equilibrium  case  (Maxwellian),  and  although  the  assumption 
of  isotropy  is  reasonable  we  do  not  know  the  actual  distribution  very  well. 
However,  it  can  be  seen  that  diffusion  is  a  function  of  the  average  energy, 
which  will  not  be  greatly  affected  by  the  shape  of  the  distribution  function. 

We  may  expect  the  greatest  change  in  the  distribution  in  the  high  energy  tail 
and  in  the  low  energy  region.  This  of  course  will  greatly  affect  the  ionization 
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and  attachment  rates,  which  depend  on  the  high  and  low  energy  electrons  and 
not  on  the  average.  In  addition  we  have  made  our  calculations  for  free  electron 
diffusion.  In  an  ionized  gas  we  also  have  positive  ions  and  except  at  the  lowest 
densities  we  do  not  expect  the  electrons  to  diffuse  freely,  but  rather  we  expect 
space  charge  forces  to  cause  the  ions  to  be  dragged  along  with  the  electrons. 
This  is  known  as  ambipolar  diffusion  and  its  rate  is  much  slower  then  free 
diffusion,  and  it  is  discussed  in  Allis  (Ref.  29), 

3)  Ionization  and  Attachment 

We  are  left  with  only  two  mechanisms  that  directly  affect  the  electron 
density,  these  are  ionization  and  attachment.  Excitation  is  important  as  a 
secondary  mechanism,  which  will  affect  the  electron  velocity  distribution. 

Since  ionization  and  attachment  are  a  function  of  the  velocity  distribution 
excitation  will  indirectly  affect  these  coefficients.  However,  we  are  using 
a  semi -empirical  approach  that  does  not  include  the  excitation  mechanism 
explicitly.  Ionization  can  occur  by  varied  mechanisms,  such  as  photo¬ 
ionization,  high  field  ionization,  and  collisions  with  excited  atoms  or  ions. 

Photo -ionization  may  be  important  in  high  pressure  streamers,  but  we  need 
not  consider  it  here.  Collisions  with  excited  atoms  (Penning  Effect)  are 
important  in  some  rare  gas  mixtures  and  may  be  important  in  a  gas  discharge 
if  impurities  are  present.  It  is  not  important  unless  excited  states  for  the 
impurity  exist  with  energy  levels  greater  then  the  ionization  potential  of  the 
main  constituent,  therefore  we  consider  only  electron-atom  collision  processes. 

Ionization  by  collision  is  essentially  a  simple  process,  by  which  all  atoms 
and  molecules  may  be  ionized.  There  is  no  difficulty  satisfying  energy  and 
momentum  considerations,  because  there  are  two  free  electrons  after  collision. 
We  can  easily  see  that  although  the  electron  collision  frequency  is  high  the 
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probability  that  an  atom  will  collide  with  an  electron  is  much  lower,  due  to 
the  low  electron  density.  We  may  therefore  consider  ionization  to  occur  in  a 
single  step,  except  for  those  atoms  that  have  a  very  long  lived  metastable 
state.  Even  then  it  is  unlikely  the  ionization  will  occur  by  multiple  collision 
since  atom-atom  collisions  will  always  be  more  probable  then  atom-electron 
collisions.  The  probability  of  ionization  upon  impact  is  zero  for  all  electron 
energies  below  the  ionization  energy.  It  rises  with  increasing  energy  until 
a  maximum  is  reached  at  3-4  times  the  ionization  energy.  Multiple  ionization 
is  possible  but  not  probable  with  the  low  ion  densities  that  exist  before  break¬ 
down.  A  great  deal  of  data  and  computation  is  available  on  ionization  potentials, 
cross -sections  and  probability,  and  Townsend  first  coefficients.  These  may 
be  found  in  the  literature  and  in  books  such  as  those  by  Brown,  (Ref.  8) 

Loeb,  (Ref.  5)  Handbuch  der  Physik  V  XXI  and  XXII,  Franklin  and  Field  (Ref.  9) 
and  in  the  recent  work  edited  by  Bates.  (Ref.  10) 

Electron  attachment  is  a  more  special  process  and  only  certain  elements 
can  form  stable  negative  ions.  To  form  a  stable  ion  the  energy  after  attach¬ 
ment  must  be  less  than  that  of  the  neutral  state.  Therefore  those  elements 
with  nearly  filled  shells  should  have  the  highest  electron  affinities  (energies 
of  attachment).  Adding  an  electron  to  an  outer  shell  is  similar  to  forming  an 
excited  state  of  the  atom,  but  the  similarity  is  not  exact  because  there  is  an 
excess  charge  which  will  interact  with  the  nucleus.  Quantum  mechanical 
calculations  have  been  attempted  for  some  of  the  simpler  atoms  for  which 
wave  functions  are  available.  (Ref.  11) 

Because  attachment  is  a  capture  process  energy  considerations  are 
important.  The  direct  capture  process 

A  +  e  - A 
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is  essentially  a  resonance  process.  Direct  capture  by  molecules  is  also 
possible  and  important.  Another  process  is  dissociative  attachment: 

AB  +  e  — >  A~  +  B 

The  process  provides  a  ready  mechanism  for  energy  and  momentum  conser¬ 
vation  and  is  likely  to  be  very  efficient.  There  is  also  the  possibility  of  a 
non-capture  process. 


AB  +  e  — A  +  B  +  e 

which  is  of  interest  in  the  formation  of  negative  ions,  but  not  as  a  loss 
mechanism  for  electrons  in  breakdown.  A  list  of  atomic  negative  ions  of 
interest  in  gas  breakdown  is  given  in  Table  I.  (Taken  from  Branscombe 
Ref.  12). 

There  are  many  molecular  attaching  gases  that  are  very  important  and 
efficient.  Buchel'nikova  (Ref.  21)  gives  the  following  information,  see 
Table  II.  The  cross  sections  in  Table  II  do  not  agree  exactly  with  certain 
other  measurements  especially  for  for  which  Craggs  (Ref.  22)  gets 
2.  3  x  10" Some  of  these  molecular  species  form  more  then  one  negative 
ion  depending  on  energy.  For  example  in  oxygen  we  see  the  following 

Oz  +  •  02- 

and 

02  +  e  +  4.  35  ev  — O  +  O 
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Table  I 

BINDING  ENERGIES  OF  NEGATIVE  IONS 


Element 

Binding  Energy 
(ev) 

Reference  and 
Measurement 

Method 

H" 

0.  8  +  0.  1 

Surface  ionization^  ^ 

He" 

0 

(14) 

Charge  exchange 

O" 

1.  465  +  0.  005 

(15) 

Photodetachment  ' 

f" 

3.  62  +  0.  09 

Surface  ionization^^ 

3.  47 

(17) 

Surface  ionization  ' 

s~ 

2.  07  +  0.  07 

(18) 

Photodetachment' 

2.  27 

Surface  ionization^ ^ 

CL‘ 

3.  76  +  0.  09 

Surface  ionization^^^ 

3.  71 

Surface  ionization^^ 

Br" 

3.  51  +  0.  06 

Surface  ionization^^^ 

3, 53  +  0.  12 

Photoionization 

r 

3.  17  +  0.  05 

Surface  ionization^**) 

3.  23 

Surface  ionization 

3.  13  +  0.  12 

Photoioni  zation  ^ 
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fl 
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fl 
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Table  II 

ATTACHMENT  CROSS  SECTIONS 


Molecule 

Cross  Section 

Energy 

cm2 

(ev) 

SF6 

5.  7  x  10"16 

0.00 

ccl4 

1.  3  x  10-16 

0.02 

cf3i 

7.  8  x  10'17 

0.05 

ccjl2f2 

-  17 

5.4  x  10  11 

0.  15 

bcl3 

- 17 

2. 8  x  10 

0.4 

HBr 

-  17 

5.  8  x  10 

0.  5 

HCL 

3. 9  x  10'18 

0.6 

h2o 

4. 8  x  10‘18 

6.4 

°2 

1.  3  x  10"18 

6.2 
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Each  of  these  has  a  different  cross  section  and  energy  dependence.  In  a 
complex  molecule  such  as  SF^  all  sorts  of  products  have  been  identified: 

(Ref,  23)  SF^",  SFg“,  SF^"  etc.  and  F~,  F^".  The  most  probable  reaction 
involves  SF^"  production,  which  has  a  sharp  peak  at  0.  03  ev.  (Ref.  24) 

SFg~  formation  has  a  much  broader  peak  at  0.  2  ev  and  a  much  lower  cross- 
section. 

It  has  been  recently  determined  (Ref.  47)  that  SF^~  is  not  completely 
stable  but  has  a  short  lifetime  (lyMsec)  after  which  it  c  ssociates  into  various 
products  such  as  SF^  +  e,  SF^”  +  F  etc.  Some  of  the  dissociation  products 
yield  free  electrons  and  some  do  not.  We  know  that  SF^  has  a  large  cross 
section  and  is  an  effective  inhibitor  of  breakdown.  Because  of  this  metastable 
character  of  SF^~  we  should  expect  different  behavior  for  short  pulse  break¬ 
down  (i.  e.  pulse  length  less  than  the  lifetime  of  the  state).  This  phenomena 
should  be  investigated  not  only  in  sulpher  hexafloride,  but  in  other  polyatomic 
gases  which  might  also  display  this  metastable  character.  In  fact  one  should 
evaluate  many  gases  with  the  consideration  of  metastable  negative  ion 
formation  since  most  previous  experiments  have  been  done  under  continuous 
or  long  pulse  conditions.  CO-,  forms  negative  ions  by  the  following  reaction 

COz  +  e  — **  CO  +  0~  , 

- 19  2 

with  a  peak  cross  section  of  5  x  10  cm  at  8  ev. 

We  see  that  only  certain  gases  will  attach  electrons,  and  that  they  do 
this  in  a  fairly  complex  manner  (i.  e.  the  functions  giving  the  energy  dependence 
of  the  cross  section  are  generally  peaked  and  not  representable  by  a  simple 
function  as  are  ionization  potentials  to  a  good  degree  of  approximation). 
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This  is  due  to  the  quantum  mechanical  nature  of  the  process.  Further  study 
is  required  on  attachment  processes  and  their  coefficients.  It  may  be  possible 
to  formulate  a  mixture  of  gases  that  has  a  high  attachment  cross  section  over 
a  wide  energy  range.  This  gas  would  then  be  an  ideal  dielectric.  In  con¬ 
sidering  breakdown  in  free  air  we  must  consider  the  presence  of  water  vapor 
which  attaches  electron  readily. 
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III.  THE  THEORY  OF  PULSE  BREAKDOWN 


We  shall  first  discuss  breakdown  that  occurs  during  a  single  pulse. 

We  assume  that  an  equilibrium,  quiesent  condition  exists  before  the  applica¬ 
tion  of  the  pulse.  We  shall  set  up  the  equations  governing  the  electron 
buildup  and  then  find  solutions  for  the  electric  field  that  will  c^tuse  the 
electron  density  to  reach  a  critical  value  in  a  time  equal  to  the  pulse  length. 

After  consideration  of  the  physical  processes  involved,  which  are 
discussed  above,  we  can  reduce  equation  (5)  to 

-J?  ■  <24  -  *a>  n  (12) 

The  source  term  S  has  been  dropped  because  for  short  pulse  breakdown  where 
the  net  ionization  is  very  large  it  will  be  important  only  as  an  initial  source 
of  electrons.  It  may  be  replaced  by  assuming  an  initial  value  for  n.  The 
spatial  dependence  has  been  removed  by  the  elimination  of  the  diffusion  term 
and  by  assuming  that  2^  and  24  are  independent  of  position.  This  is 
equivalent  to  assuming  that  the  electron  distribution  function  is  isotropic 
and  homogeneous  in  space.  The  ionization  and  attachment  rates  are  average 
values  over  the  electron  distribution.  The  solution  to  equation  (12)  is 

Z* 

n  =  nQ  e  , 

where 

V'Vi-K  ■ 
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These  concepts  and  general  approach  have  worked  for  c.  w.  and  long 
pulse  breakdown  (Refs.  35,  37,38)  where  diffusion  losses  were  included.  If 
the  neglect  of  diffusion  is  justified,  and  this  approximation  is  better  for 
shorter  pulses,  this  approach  should  be  valid  for  short  pulse  breakdown. 

Only  if  there  were  processes  which  are  important  for  very  short  times  would 
we  expect  a  discrepancy.  Such  effects  as  metastable  states  of  negative  ions 
are  not  expected  to  affect  the  theory,  they  just  require  re-evaluation  of  the 
constants. 

Let  us  examine  the  theory  for  possible  limitations  in  the  very  short 
pulse  limit.  By  assuming  constant  collision  frequency  we  see  that  there  is  a 
lower  limit  on  the  time  between  collisions,  and  that  if  the  pulse  length 
is  less  than  this  time  there  can  be  no  ionizations.  This  may  be  expressed 
as^^c  ^  **  Using  our  value  of  collision  frequency  for  air  reduces  this 
to  pT<*0.  2  nanosecond  torr  for  collisions,  and  hence  ionizations  to  occur. 

This  puts  a  lower  limit  on  p  of  0.01  torr  for  a  20  nsec  pulse.  This  theory 
is  not  exact  because  the  collision  frequency  is  a  function  of  energy  and  there 
will  be  a  distribution  of  collision  times.  By  considering  somewhat  more 
realistic  potentials  for  interaction  of  electrons  and  atoms  (Ref.  26)  one  may 
see  that  may  decrease  with  increasing  energy.  This  can  lead  to  run-away 
electrons  because  as  the  collision  frequency  decreases  we  must  increase  the 
energy  ot  insure  a  sufficient  number  of  ionizations.  This,  however,  reduces 
the  number  of  collisions,  and  we  see  that  breakdown  will  not  occur.  This 
effect  will  be  very  sensitive  to  the  interaction  potential. 

As  the  mean  free  path  increases  with  decreasing  pressure  more  and 
more  electrons  will  be  swept  to  the  walls  at  each  cycle.  This  causes 
secondary  effects  which  are  discussed  by  Brown.  (Ref.  4)  For  the  higher 
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pressures  and  this  is  not  important.  Therefore  as  long  as  pt^O.  2 

and  p  is  not  too  low  we  should  expect  the  theory  to  hold. 

We  must  now  investigate  the  calculation  of  To  do  this  in  detail 
from  elementary  principles  requires  a  knowledge  of  the  distribution  of  electron 
energies,  the  attachment  and  ionization  cross  sections  as  a  function  of  energy, 
and  the  collision  frequency  as  a  function  of  energy.  Even  if  we  had  all  this 
information,  we  would  still  have  an  extremely  difficult  problem.  To  find  the 
distribution  function  one  must  solve  the  Boltzmann  equation 

4j4  +\7  *  Vrf+  £  ‘  Vvf  =  (-JT>  coll  <13) 

The  collision  term  on  the  right  involves  the  effects  of  all  collisions  (elastic, 
exciting,  attaching  and  ionizing)  which  in  general  not  known  exactly  and 
usually  varies  in  a  complex  manner.  Brown  and  Allis,  (Ref.  27)  and  Desloge 
and  Matthysse  (Ref.  28)  have  given  approximations  for  this  term.  A  general 
solution  to  this  is  virtually  impossible  and  in  electrical  discharge  problems 
one  usually  expands  f  in  spherical  harmonics  (see  Allis  (Ref.  29)  or 
Delcroix  (Ref.  26)).  Even  then,  however,  a  solution  is  not  easily  obtained 
except  by  approximate  methods  and  usually  only  with  great  difficulty.  The 
Boltzmann  equation  approach  has  been  used  to  give  a  theory  of  breakdown 
for  hydrogen,  (Ref.  30)  helium,  (Ref.  31)  helium  without  excited  states,  (Ref.  32) 
and  neon.  (Ref.  33)  To  apply  this  to  other  gases,  especially  to  the  more 
complex  molecular  species,  would  be  very  difficult.  In  any  case,  there  is  an 
extreme  lack  of  cross-section  information  for  many  of  the  gases  in  which  we 
have  interest.  The  Boltzmann  equation  may  also  be  used  in  a  more  general 
way  to  establish  the  importance  of  relaxation  times  for  the  various  collision 
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processes  and  to  gain  insight  into  the  effects  of  certain  perturbations. 

The  usual  study  of  breakdown  or  ionization  processes  avoids  calculation 
from  first  principles  although  some  work  of  this  kind  has  been  done.  (Ref.  34) 
The  theory  usually  used  has  been  semi-empirical  utilizing  empirically 
determined  ionization  coefficients  and  theoretically  supplying  the  other  details 
of  the  calculation.  This  is  a  very  useful  method  and  further  work  along  these 
lines  is  suggested.  Empirical  determination  of  the  constants  has  the  advantage 
that  one  automatically  obtains  averages  over  the  distribution  and  that  the 
difficult  problem  of  accounting  for  excitation  is  automatically  included  in  the 
measurement.  The  only  problem  then  is  to  determine  whether  the  long  time 
c.  w.  measurements  of  ionization,  and  attachment  differ  from  very  short 
pulse  measurements  because  of  the  finite  time  it  takes  to  establish  the  steady 
state  distribution.  This  also  depends  on  the  electron -atom,  and  electron- 
electron  interactions  and  varies  with  the  form  of  the  potentials.  Delcroix 
(Ref.  26)  has  a  discussion  of  the  problem,  and  calculates  relaxation  times  for 
two  cases.  Except  for  very  small  values  of  p  X  we  expect  that  this  effect 
will  not  be  important. 

The  way  to  evaluate  breakdown,  then  is  to  use  experimental  data. 

For  air  Gould  and  Roberts  (Ref.  35)  used  d.  c.  coefficients  and  an  equivalent 
field  method  to  determine  microwave  coefficients.  This  approach  may  lead 
to  errors,  because  there  is  an  a  priori  assumption  made  that  the  distribution 
functions  are  the  same  for  the  equivalent  values  of  E/p.  This  is  not  known 
to  be  true  and  the  agreement  may  be  fortuitous.  Holstein  (Ref.  36)  has, 
however,  shown  that  electron  distributions  in  microwave  discharges  approach 
those  in  the  d.  c.  case  for  a  limited  range.  The  method  of  Gould  and  Roberts 
may  be  used  to  calculate  breakdown  for  gases  in  which  microwave  data  is  not 
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available,  but  for  which  Townsend  measurements  have  been  made.  Two 
other  papers  on  breakdown  air  by  MacDonald,  (Ref.  37)  and  by  Kelly  and 
Margenau  (Ref.  38)  both  utilize  Brown's  (Ref.  39)  experimental  microwave 
ionization  data  obtained  from  breakdown  measurements.  This  coefficient 
is  a  microwave  coefficient  for  the  diffusion  limited  case  in 
analogy  to  the  mobility  limited  Townsend  ionization.  The  approach  in  these 
papers  is  somewhat  different,  but  the  essential  results  are  the  same.  As 
before,  we  have  the  equation 

4?  •  <H  -  va) »  ,  <«> 

and 

n(t)  =  no  2* 

We  define  breakdown  such  that  the  electron  density  at  the  end  of  a  pulse  is 
to  be 

nb(£)  =  1.24  x  10‘8  f2  (15) 

so  that  breakdown  will  occur  if 

>*=4.6  log  f -26.0  (16) 

3  3 

where  n  was  chosen  as  3  x  10  /cm  • 
o 
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Kelly  and  Margenau  utilize  the  modified  Boltzmann  equation  and 


introduce  a  distribution  function 


J(v,  t)  =  0(v) 


(17) 


The  only  result  of  this  is  to  formally  show  that  D  =  .which  was 

somewhat  implicit  in  their  choice  of  the  distribution  function. 

We  have  a  combination  of  theory  and  empirical  data  which  will  allow 
us  to  calculate  the  breakdown  fields.  We  can  calculate  the  value  of 
required  to  cause  breakdown,  as  a  function  of  frequency,  with  equation  (16). 

a  2 

Herlin  and  Brown  (Ref.  39)  have  tabulated  y  as  a  function  of  (E/p )  . 
Therefore  if  one  knows  D  as  a  function  of  E/p  it  is  possible  to  relate  the  net 
ionization  "JJ  and  the  field  strength  E.  Unfortunately  T  is  also  a  function  of 
p  \  (Xio  the  wavelength  which  corresponds  to  the  frequency-^  ),  however 
when  pj^  is  greater  than  200  all  the  curves  merge.  For  p  \  200  it  is 
convenient  to  write 


Vt  =  Dp  *  (|-)2>  pf.  .  <18> 

which  is  obtained  from  Herlin' s  coefficient  by  judiciously  multiplying  by  the 
correct  factors.  For  smaller  values  of  pA  we  must  multiply  equation  (18) 
hy  A  »  but  for  the  higher  pressures  we  can  obtain  a  universal  breakdown 
curve  from  equation  (18).  All  the  factors  are  either  known  or  may  be  calcu¬ 
lated.  The  only  difficulty  is  with  Dp.  We  have  previously  shown  that 

Dp  =  2  x  10^  ^u^  . 
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The  difficulty  is  in  evaluating  .  Margenau  uses  a  distribution  function 

that  he  has  previously  (Ref.  40)  shown  to  be  approximately  correct,  whereas 
MacDonald  uses  the  (u^  data  of  Healy  and  Read,  (Ref.  41)  which  was  taken 

in  dc  discharges.  Margenau' s  distribution  is  practically  independent  of 

6  4 

E/p,  and  gives  Dp  =  1.6  x  10  ,  whereas  MacDonald  obtains  Dp  &  1.  2  x  10  E/p 

(for  higher  pressures).  These  will  be  equal  for  E/p  =  130,  wtych  means  that 
their  calculations  of  Dp  are  about  equal  in  the  range  where  breakdown  occurs. 
This  is  of  course  to  be  expected,  because  the  data  for  calculation  of  these 
values  is  supposedly  based  on  actual  measurements  in  gas  breakdown  or  at 
least  on  a  theory  that  is  known  to  yield  proper  values.  The  agreement  is, 
however,  far  from  perfect  and  we  expect  that  this  discrepancy  may  be  over¬ 
come  by  attempting  to  measure  diffusion  losses  directly.  Although  the  diffusion 
coefficient  is  not  sensitive  to  the  shape  of  the  distribution  it  does  depend 
directly  on  the  value  of  the  energy,  which  in  turn  must  be  known  absolutely. 
Unfortunately  we  do  not  generally  know  the  average  energy  either  directly  or 
indirectly. 

Using  equation  (18)  and  a  value  of  Dp  =  1.  36  x  10^  we  have  plotted  a 
universal  breakdown  curve  (Figure  2).  This  curve  agrees  with  previous 
calculations  of  Gould  and  Roberts  (Ref.  35)  and  at  the  longer  pulse  lengths 
agrees  with  the  published  measurements  of  MacDonald.  This  curve  can  be 
used  to  predict  breakdown  in  planning  experiments,  so  that  one  may  predict 
the  requirements  of  the  equipment  needed  to  do  pulse  measurements. 

2)  Multiple  Pulse  Breakdown 

It  is  possible  over  a  series  of  pulses  to  build  up  the  electron  density 
so  that  breakdown  occurs  after  N  pulses.  We  can  write  the  following  two 
equations  for  this  situation, 
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The  buildup  equation  during  each  pulse  is 

-gy  =  *Pn  *  (12) 

and  the  decay  equation  between  pulses  is 

-^r  =  •  34  n  "  <Kn2  +  D72  n  (18) 

We  include  the  recombination  and  diffusion  in  the  decay  because  of  the  longer 
time  between  pulses.  In  a  strongly  attaching  gas  we  may  simplify  equation  (18) 
to 


(18a) 


The  electron  density  after  N  pulses  is 

N*ti  -  <N  -  1)24 tv 

n(t)  =  nQ  e  (19) 

or 

lnJL=  -  (N  -  l)12at-2  (19a) 

o 

There  is  a  -  term  in  but  it  will  have  a  different  value  than  during 
the  decay  period  because  the  electrons  in  the  long  period  between  pulses 
will  come  into  equilibrium  with  the  gas  and  hence  will  have  a  Maxwellian 
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distribution  at  a  temperature  of  about  300°K.  The  attachment  cross  section, 
is  high  for  low  energy  electrons  so  that  the  attachment  during  decay  will*  bn 
an  efficient  process.  Unfortunately  in  the  microwave  measurements  of  ^ 
there  is  no  way  to  separate  the  ionization  and  attachment  terms.  Either  d.  c. 
data  will  have  to  be  used,  or  a  direct  measurement  may  have  to  be  made  in 
the  afterglow.  This  should  be  a  prime  consideration  in  the  determination  of 
new  coefficients. 
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IV.  MICROWAVE  TESTING  METHODS 

Experimentally  we  require  a  faet  rising  pulse  that  can  supply  the 
required  field  strength.  The  field  strength  is  related  to  power  density  so  it 
is  a  matter  of  getting  fast  pulses  and  power.  We  do  not,  for  purposes  of 
testing,  need  actual  short  pulses  although  to  supply  the  power  we  need  in  a 
long  pulse  will  be  difficult.  We  may  divide  this  problem  into  two:  the  problem 
of  rise  time,  and  the  problem  of  field  strength. 

We  may  relate  the  required  field  strength  to  power  density  by  the 
Poynting  theorem 


S  =  E  x  H  (20) 

By  specifying  the  geometry  we  may  relate  E  and  H  through  an  impedance  for 
simple  modes.  The  following  cases  are  particularly  useful: 

The  Plane  Wave  in  Free  Space 

p  -  1  E2 

P  “7  377 

Rectangular  Wave  Guide 

TEjq  mode 

where  a  and  b  are  the  dimensions  of  the  guide  and  Z  is  the  impedance 
(about  500  ohms). 

We  cannot  use  a  normal  section  of  waveguide  as  a  test  section  for  two 
reasons.  First,  we  must  handle  the  same  power  everywhere  so  that  we  can 
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have  a  breakdown  anywhere  in  the  system,  and  second,  because  the  guide  is 
parallel  we  cannot  predict  where  the  breakdown  would  occur.  Furthermore, 
to  reduce  our  peak  power  requirement  we  wish  to  intensify  the  fields.  The 
fields  in  a  resonant  cavity  are  of  course  very  high,  and  this  is  used  in  c.  w. 
breakdown  for  the  evaluation  of  the  ionization  coefficients.  However  high  Q 
cavities  have  a  long  time  constant,  t  =  2  Q/(jJ,  which  is  of  the  order  of  a 
microsecond  for  X  band  cavities.  This  makes  the  cavity  unsuited  for  pulsed 
measurement. 

A  simple  method  for  obtaining  higher  fields  is  to  taper  the  waveguide. 

In  the  TEjq  mode  the  propagation  is  not  affected  by  the  height*  so  that  we 

can  reduce  the  height  and  increase  the  power  density  and  hence  the  field 

strength.  The  field  strength  at  constant  power  is  inversely  proportional  to 
2 

b  .  A  further  gain  may  be  realized  by  using  ridged  waveguides.  The  complete 
equations  for  power  and  field  distribution  (Ref.  42)  are  somewhat  involved, 
and  one  may  examine  the  original  papers  for  a  complete  description*  However, 
we  may  calculate  the  field  strength  in  the  following  approximate  but  simpler 
manner.  Consider  a  ridge  waveguide  with  the  following  dimensions  (Figure  3). 
From  Hopfer  (Ref.  43)  we  find  that  ’Xc  =  3.  5  A  compared  to  =  2  a  for 
rectangular  guide,  therefore  we  may  reduce  A  to  0.  8a,  and  S  will  be  0.  4a, 

B  will  be  0.  8b  and  D  will  be  0.  24b.  The  fields  are  concentrated  in  the  region 
of  the  ridge  so  that  we  see  that  much  greater  power  densities  are  experienced. 
Hopfer  gives  information  that  also  allows  one  to  calculate  the  field.  For  this 
case  we  find  that  the  guide  would  have  E/p  =  40  with  28  KW,  whereas  a  tapered 
rectangular  guide  with  the  same  actual  height  would  require  108  KW,  This 
is  only  an  approximate  calculation  but  shows  the  advantage  in  the  ridged  guide. 
The  ridge  guide  may  be  tapered.  To  provide  a  transition  from  rectangular 
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guide  it  must  be  tapered  in  both  dimensions  (that  is  the  width  and  the  ridge 
gap).  There  is  a  minor  problem  with  fringing  fields  at  the  corners,  but  this 
can  be  eliminated  by  contouring  the  ridge.  The  tapered  ridge  guide  provides 
a  convenient  mechanism  for  providing  high  fields  with  modest  power  require¬ 
ments;  it  also  provides  a  localized  breakdown  because  the  rectangular  guide 
that  feeds  the  ridged  section  will  have  much  lower  fields.  Loss  of  electrons 
by  diffusion  may  have  to  be  included  for  longer  pulses  because  the  diffusion 
length  is  small. 

There  is  also  the  possibility  of  using  mismatched  sections  to  produce 
higher  fields.  If  the  test  section  were  placed  a  quarter  wavelength  in  front 
of  a  short  we  expect  the  following  sequence.  A  traveling  wave  propagates 
down  the  guide  toward  the  test  point.  The  amplitude  builds  up  in  a  time  equal 
to  the  rise  time,  but  there  will  be  a  backward  wave  reflected  from  the  short 
which  will  set  up  a  standing  wave  giving  a  E  field  maximum  at  the  test  point. 
Since  the  propagation  time  is  very  short  (over  a  half  wavelength  path)  in 
comparison  to  the  rise  time,  we  expect  the  field  to  rise  to  2  Eq,  with  a  time 
characteristic  of  the  pulse  rise  time.  The  time  between  the  arrival  of  the 
traveling  wave  and  the  establishment  of  a  standing  wave  is  extremely  short 
CV10-10  sec).  We  must  remember  that  in  a  wave  guide  the  pulse  shape  will 
be  degraded  due  to  dispersion  if  long  propagation  paths  are  used  so  that  we 
must  keep  our  wave  guide  length  at  a  minimum. 

We  also  considered  the  production  of  fast  rising  pulses.  It  appears 
that  only  an  experimental  program  can  provide  new  solutions  to  this  problem. 
Presently  available  magnetrons  have  pulse  rise  times  in  the  order  of  20-50  nsec; 
limited  by  mode  changing  when  faster  drive  pulses  are  used.  It  has  been 
reported  (Ref.  48)  that  special  designs  to  eliminate  mode  changing  can  lower 
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the  rise  time  to  say  1.  0  nsec.  This  has  not  been  fully  exploited  because 
present  rise  times  are  probably  satisfactory  with  microsecond  pulses.  We 
can  of  course  use  a  fast  rise  long  pulse  for  testing  if  a  source  were  available. 
It  is  unlikely  that  we  can  achieve  the  high  peak  powers  required  for  breakdown 
in  nanosecond  times  in  a  long  pulse. 

In  order  to  get  fast  rise  time  one  can  consider  either  pulse  shaping  or 
a  new  generating  scheme.  If  the  wave  guide  could  be  suddenly  shorted  while 
a  wave  were  propagating  in  it  we  would  expect  a  sharp  reflection.  A  useful 
fast  switch  is  a  gas  discharge  (e.  g.  ,  the  TR  Tube).  We  have  seen  that  for 
U)y"^c  *ke  reflection  as  the  electron  density  approaches  the  plasma  resonant 
density  is  rather  abrupt,  so  that  we  expect  a  sharp  reflection  from  a  rapidly 
ignited  low  pressure  discharge.  This  could  provide  the  means  for  pulse 
shaping.  A  circulator  could  be  used  to  distribute  the  pulse  to  the  test  section 
and  keep  the  reflected  and  incident  waves  separated.  This  system  appears 
feasible  and  depends  mainly  on  the  fast  switching  which  is  available.  We  have 
been  informed  by  W.  Quinn  of  RADC,  however,  that  this  has  been  tried  and 
has  not  been  very  satisfactory.  This  kind  of  pulse  shaper  would  depend  on 
the  fast  switch,  the  magnetron  pulse  and  peak  power  available  and  on  the 
properties  of  the  circulator.  It  might  be  possible  to  develop  this  system,  but 
there  has  appeared  in  the  literature  (Ref.  45)  recently  a  more  promising 
device,  the  variable -coupling  ring  resonator  which  can  provide  fast  rise, 
short,  high  peak -power  pulses.  This  device  consists  of  a  resonant  re-entrant 
wave  guide  coupled  to  a  straight  section  of  guide.  (Figure  5).  This  device 
has  been  previously  used  as  a  field  intensity  magnifier,  and  as  an  elimination 
filter.  An  electric  field  intensity  gain  in  the  ring  of  about  10  can  be  obtained, 
but  because  it  is  a  high  Q  device  the  rise  time  is  too  long  to  make  it  useful 
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to  us  as  a  method  of  obtaining  high  fields.  However,  if  after  the  ring  is 
charged  the  coupling  coefficient  of  the  ring  is  suddenly  changed  to  unity,  it  is 
possible  to  take  out  all  the  energy  in  a  single  pulse.  The  rise  time  is  limited 
by  the  switch  in  the  coupler  and  the  pulse  length  depends  on  the  ring  length. 
Rise  times  of  3  nsec  with  pulse  length  of  20  nsec  and  a  peak  power  of  750  KW 
have  been  achieved.  (Ref.  45)  The  switching  again  depends  on  a  rapidly  ignited 
gas  discharge.  This  appears  to  be  thj  best  method  presently  available  for 
obtaining  short,  fast  high  energy  pulses,  and  would  certainly  make  a  useful 
device  for  breakdown  experiments. 
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V.  SUGGESTIONS  FOR  FUTURE  WORK 


One  of  the  main  aims  of  this  project  was  the  investigation  of  the 
desirability  and  the  requirements  for  doing  pulse  breakdown  experiments. 

We,  have  shown  in  proceeding  sections  that  a  theory  exists  for  predicting 
breakdown  fields.  This  theory  is  essentially  phenomenological  and  requires 
the  use  of  empirically  determined  coefficients.  This  is  entirely  suitable  for 
the  guidance  of  system  designers,  but  it  does  indicate  a  lack  of  fundamental 
knowledge  of  the  processes  involved  in  microwave  breakdown.  It  is  extremely 
desirable  to  have  a  detailed  microscopic  understanding  of  the  processes 
involved  in  breakdown,  and  ultimately  this  approach  is  best.  However,  there 
are  inherent  difficulties  in  the  microscopic  approach,  and  we,  therefore 
suggest  the  following  three  part  program,  which  is  designed  to  give  us  the 
data  on  breakdown  that  is  required  for  system  design,  and  to  give  a  more 
basic  understanding  of  the  processes. 

1)  Pulse  Breakdown  Measurements 

A  primary  requirement  is  for  actual  measurements  under  pulsed 
conditions.  We  suggest  a  series  of  pulse  breakdown  measurements  in  at 
least  the  following  gases:  air,  humid  air,  CC>2>  SF^,  and  a  none  attaching 
rare  gas.  These  measurements  would  be  made  using  presently  available 
techniques  such  as  the  variable -couplding  ring  resonator,  along  with  tapered 
wave  guide.  The  suggested  gases  were  chosen  because  they  represent  the 
important  cases  encountered  in  present  systems  and  because  they  represent 
a  range  of  attachment  cross  sections  (see  Table  II).  The  rare  gas  or  other 
non-attaching  gas  is  suggested  as  a  check  on  the  approximations  made  in  the 
theory.  If  diffusion  and  recombination  are  negligible  and  attachment  is 
absent  then  there  are  no  loss  mechanisms  and  breakdown  should  be  controlled 

ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 

31  ARF-A2 17-7 

— . *s>  ■  ~  Ll:  ■> 


entirely  by  the  ionization  coefficient.  This  direct  measurement  of  the  break¬ 
down  fields  is  not  enough  since  it  gives  only  data,  and  does  not  provide  an 
explanation  of  the  phenomena.  To  provide  the  understanding  we  must  perform 
supplementary  experiments. 

2)  Measurement  of  Ionization  and  Attachment  Coefficients 

We  have  a  theory  of  microwave -pulsed  breakdown  that  requires  only 
that  we  insert  the  proper  coefficients.  The  ionization  coefficient  can  be 
obtained  from  c.  w.  breakdown  measurements  as  originally  proposed  by 
Herlin  and  Brown.  (Ref.  6)  This  involves  the  study  of  breakdown  in  a  cavity, 
and  yields  a  net  ionization  coefficient.  It  might  be  possible,  by  using  a  second 
signal  as  a  probe,  to  determine  electron  density  for  fields  less  then  that 
required  for  breakdown.  This  in  turn  would  yield  the  net  ionization  coefficient 
for  low  values  of  E/p,  and  we  therefore  might  be  able  to  separate  the  attach¬ 
ment  and  ionization  coefficients  in  a  manner  similar  to  that  used  in  the  analysis 
of  d.  c.  discharges. 

A  similar  cavity  measurement  can  also  be  used  to  measure  the  net 
decay  coefficient  in  the  afterglow.  In  this  case  it  is  possible,  from  the  shape 
of  the  curves,  to  separate  (at  least  partially)  the  attachment  and  recombination 
coefficients.  These  two  coefficients,  net  growth  and  net  decay,  would  allow 
us  to  predict  breakdown  fields  for  widely  varying  conditions  under  both  single 
pulse  and  multi  pulse  conditions.  Comparison  of  such  predicted  results  with 
the  direct  pulse  breakdown  measurements  would  test  the  accuracy  of  the 
theory.  If,  as  we  believe,  the  theory  is  sufficiently  accurate,  then  we  could 
substitute  c.  w.  measurements  for  pulse  measurements  and  still  predict 
breakdown.  This  would  be  extremely  advantageous  since  it  would  eliminate 
a  large  number  of  tedious  measurements  over  a  wide  range  of  parameters 
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such  as  pulse  width. 

This  still  does  not  contribute  to  our  understanding  of  the  microscopic 
processes*  in  order  to  do  this  we  must  examine  a  more  detailed  theory. 

3)  Microscopic  Theory  of  Breakdown 

By  a  microscopic  theory  we  imply  one  in  which  the  details  of  the 
collision  processes  are  used.  This  requires  a  knowledge  of  the  electron 
velocity  distribution  and  detailed  knowledge  of  the  cross  sections  for  the 
various  processes.  If  one  had  this  detailed  information  it  would  be  possible 
to  calculate  the  coefficients  discussed  in  the  preceeding  section.  With  a 
somewhat  less  detailed  knowledge  it  is  possible  to  compare*  or  calculate, 
the  microwave  ionization  from  the  Townsend  first  coefficient. 

What  we  propose  therefore  is  to  attempt  solutions  to  the  Boltzmann 
equation  in  order  to  improve  the  approximation  to  the  distribution  function. 
This  will  be  used  to  calculate  ionization  or  attachment  coefficients  for  cases 
where  the  cross  sections  are  known.  These  calculations  will  be  compared 
to  the  microwave  measurements.  The  aim  being  to  use  both  the  theoretical 
calculations. 

If  one  has  available  a  function  of  v,  g(v)  and  its  average  Jg(v)  jn(  v)  dv, 
it  is  not  possible  to  establish  the  distribution  function  £(v),  but  it  might  be 
possible  to  use  such  information  to  decide  between  various  possible  distri¬ 
butions  or  to  adjust  constants  in  a  calculated  distribution. 

The  suggested  program  realizes  the  usefulness  of  the  phenomenological 
approach  and  also  realizes  the  importance  of  providing  immediate  data  on 
breakdown.  It  further  realizes  the  importance  of  gaining  a  complete  and 
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detailed  understanding  of  the  phenomena.  In  our  suggested  program  we  have 
tried  to  stay  with  presently  available  techniques  and  concentrate  on  the 
phenomenon  of  breakdown. 
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Figure  4.  Schematic  Design  of  Pulse  Generator  Using  Circulator. 
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D  diffusion  coefficient 

e  electron  charge 

E  field  strength  (volts  /meter) 

f  frequency  (cycles/sec) 

fp  resonant  plasma  frequency 

£  velocity  distribution  function 

mean  free  path  (of  electrons) 
m  mass  of  electron 

n^  number  density  of  electrons  at  breakdown 

P  pressure  (Torr) 

R  reflection  coefficient  for  EM  waves 


S  generalized  source  term  for  electron  production 

2 

S  Poyntings  vector  (watts /m  ) 

u  electron  energy 

electron  energy  averaged  over  velocity  distribution  function 
v  electron  velocity 

(v)  electron  velocity  averaged  over  the  distribution  function 

* 

oC  recombination  coefficient  cm  /ion  sec 

P  diffusion  particle  current  density 


gradient  in  coordinate  space 
^  gradient  in  velocity  space 
yj  generalised  index  of  refraction 
^  microwave  ionization  coefficient 
^  wavelength 
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APPENDIX  I 


DETAILS  OF  THE  CALCULATION  OF  THE 
DIFFUSION  COEFFICIENT 

The  "particle  current"  density  (excluding  charge)  may  be  written  as 

r  =  -V<Dn>  •  (I- 1) 

This  equation  includes  all  the  diffusion  effects  due  to,  density  gradients, 
temperature  gradients,  and  any  other  effect  that  might  effect  D.  However 
we  are  interested  only  in  density  effects  and  may  rewrite  this  as 

P  =  -  D\/n  (I- la) 

This  can  be  thought  of  as  a  defining  statement  for  O.  In  order  to  get  a 
mathematical  expression  and  physical  interpretation  for  D  we  remember  that 
p  may  also  be  written  as 

f1  *»(v)  .  (1-2) 

For  simplicity  consider  the  one -dimensional  case.  Let  A  x  equal  the 
average  distance  a  particle  moves  between  collisions,  and  £*  the  time 
between  collisions.  Then  in  a  small  region  of  space  (x  -  A  x/2,  x  +  A  x/2) 
the  net  current  is 

[*♦**«  !$-<»-¥-!!>]•  «-»> 
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which  reduces  to 


n  -  (Ax)2  <5n 
1  x  "  £/  "gx 


(I -3a) 


Therefore  we  see  that 


D  = 


(Ax)2 

t 


d-4) 


which  shows  that  D  is  the  mean -squared-distance  a  particle  moves  in  unit 
time.  The  collision  frequency  (  2£)  is  just  equal  to  the  1/p  ,  so  that  we  may 
also  write 


d  =  (Ax)2l;  c 


(I -4a) 


A  more  exact  calculation  may  be  made  by  starting  from  the  Maxwell* Bolt smann 
equation  for  the  electrons. 


if 

U 


(1-5) 


The  implications  and  approximations  are  described  in  many  standard  works 
(c.  f.  Rose  (46)).  If  we  assume  further  that  the  distribution  function  £ (v,  r) 
is  equal  to 


J*  (V,  r)  =  (r,  v)  +  Jj  (r,  v) 


(1-6) 
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where  (r,  v)  i;$  isotropic,  and 


n(r)  =  J'dv  (r,  v)  (1-7) 

is  the  number  of  electrons.  We  assume  j  is  a  small  perturbation  on  the 
i30tropic  distribution  and  that  its  integral  over  all  phase  space  is  zero.  We 
introduce  this  distribution  function  into  equation  (1-5)  and  look  for  an 
equilibrium  solution,  whence 

V7r  <fo  +  fl>  +  K-Vv  '  fo  +  fl>  *  +  4^*  a'8) 

The  higher  order  terms  Vr  f  1*  Vv'fo  are  neglected  and  ("  ^  equals 

zero  since  it  is  by  definition  an  equilibrium  distribution.  The  (1-8)  reduces 


>*rfo4 


The  anisotropy  -j  j  is  generated  by  the  density  gradient  and  the  external  force 
F.  We  assume  that  the  distribution  function  would  relax  to  £o>  in  the  absence 
of  forcing  terms,  due  to  momentum  transferring  collisions  and  then 


(i- 10) 


where  '))  is  the  momentum  transferring  collision  frequency  (Ref.  29).  We 
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have 


P  =  n  (  =|  d  v  vf  =  d?  *  ( £0  +  £\) 

S  jdv  v  f  j  , 


(1-11) 


because  vjp  Q  is  zero  due  to  isotropy  of pQ.  We  may  now  substitute 

for  in  equation  (1-9)  dropping  the  force  term  as  we  are  interested  only 


in  diffusion  due  to  density  gradients.  Then 


r  =.  jdv  v(^-»  Vr  fQ>  •  (i-12) 

Consider  a  gradient  only  in  the  x  direction.  Then  v*  Vxf  o  ■  Vi 
and  the  flow  will  also  be  in  the  x  direction.  Therefore 


Jx_  ^fo' 


(1-13) 


2  2  2  2 
V  +  V  +  V  =  v  , 

x  y  z 


so  that  near  equilibrium  where  the  distribution  is  almost  isotropic  the 
component  velocities  are  equal  and  for  any  component  we  have 


2  m 

V 


(1-14) 
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